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We present a new approach to study the magnetic phase stability of magnetic overlayers on 
nonmagnetic substrates. The exchange integrals among magnetic atoms in the overlayer are esti- 
mated in the framework of the adiabatic approximation and used to construct the effective classical 
two-dimensional Heisenberg Hamiltonian. Its stability is then studied with respect to a large num- 
ber of collinear and non-collinear magnetic arrangements which include, as special cases, not only 
ferromagnetic and various antiferromagnetic configurations, but also possible incommensurate spin- 
spiral structures. This allows us to investigate a broader class of systems than a conventional total 
energy search based on few, subjectively chosen configurations. As a case study we consider the 
Fe-monolayer on the random nonmagnetic bcc-Ta^Wi-^OOl) surface which was studied recently 
by a conventional approach. We have found a crossover of the ground state of the Fe monolayer 
from the ferromagnet on the Ta surface to the c(2 x 2) antiferromagnet on the W surface and that 
at the composition with about 20 % of Ta an incommensurate magnetic configuration might exist. 

PACS numbers: 71.15.Mb, 75.30.Et, 75.70.Ak 

I. INTRODUCTION 

A deeper understanding of the magnetic ground state as well as of finite-temperature properties of a system with 
local magnetic moments can be obtained in terms of corresponding exchange interactions. This is particularly true 
for new, artificially prepared systems, like, e.g., monolayers of magnetic atoms on nonmagnetic substrates. A strong 
hybridization of the electronic states of the magnetic atoms and of the substrate may significantly influence the 
magnetic ground state of these systems. We mention as a typical example the Fc-monolayer on bcc-W(OOl) substrate. 
While bcc-Fe in bulk state is a ferromagnet (FM), the Fe-monolayer on the bcc-W(OOl) surprisingly has a c(2 x 2) 
antiferromagnetic (AFM) ground state, as was recently found by spin-polarized STM measurements^. 

In a subsequent theoretical study, Ferriani et al£ demonstrated a strong dependence of the magnetic state of an 
Fe-monolayer on the substrate by studying its stability on the (OOl)-facc of bcc-Ta x Wi_ x random substrates. It was 
found that the ground state of the Fe-monolayer on Ta(001) is the FM rather than the AFM one and there thus 
exists an interesting crossover between these two magnetic configurations at an intermediate substrate composition. 
Results of this study also confirm the robustness of such behavior - even the unrelaxed system exhibits the same 
behavior although large inward relaxation of about 18 % exists between the Fe monolayer and the first substrate 
W-layer, which can modify the result quantitatively. The authors of Ref. [U also clearly demonstrated the usefulness 
of exchange interactions for a deeper understanding of the phenomenon. 

Based on these facts we wish to present here a more general approach to the study of the magnetic stability of such 
systems. Let us first remind the reader of the approach adopted in paper— as well as in a number of other studies^— 
Exchange interactions are estimated by evaluating the total energies of the FM and a few AFM configurations by 
conventional bandstructure methods. In particular, in the case of Fe overlayer on the W(001) or Ta(001), the authors of 
Ref. chose c(2 x 2)- and p(2 x 1)-AFM confi gurations and mapped the total energy differences of these configurations 
and of the reference FM state onto the classical two-dimensional Heisenberg Hamiltonian (HH). In this way they 
obtained two nearest-neighbor exchange couplings. This approach has certain limitations: (i) the choice of AFM 
configurations is dictated by requirements of simplicity, although any other configuration can be used as well. One 
has to be aware that if some atypical configuration is chosen, the results can be biased; (ii) the number of determined 
exchange interactions depends on the number of chosen configurations, which is limited. This is in a striking contrast 
with the existence of long-range interactions in such systems. For example, the presence of the surface state, e.g., on 
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the fcc(lll) surfaces of noble metals leads to the decay of corresponding exchange interactions with the distance d as 
d~ 2 (see Ref. H^) as compared to the <i _3 -decay known from bulk systems (the RKKY-interactions). In addition, the 
bandstructure methods do not yield directly individual exchange interactions, but rather their certain infinite sums; 
(iii) the use of conventional bandstructure methods is strongly complicated in the presence of disorder, either in the 
magnetic overlayer or in the substrate. On the other hand, if the substrate disorder is weak, a simple approach (the 
virtual-crystal approximation) can be successfully used as it was done in Ref. [H for the case of Fe-overlayer on the 
TaW-alloy substrate. 

An alternative approach is to estimate magnetic interactions among magnetic atoms directly by evaluating the 
energy cost related to small rotations of spins at arbitarily chosen two sites with the help of the adiabatic theorem 
and the Green function method. This approach, first introduced by Lichtenstein et alJ£ for bulk ferromagnets and 
by Oguchi et al^ for bulk paramagnets within the disordered local moment (DLM) state, was then successfully gen- 
eralized to a number of magnetic systems including magnetic overlayers on crystal substrates^ (for further magnetic 
systems see the review papei— ). In this method the limitations mentioned above are removed. As the method is 
formulated in the real space, the systems without perfect translational symmetry (e.g. random alloys, overlayers, and 
surfaces) can be treated on an equal footing with the crystals. 

One remark concerns the choice of the reference state from which the exchange interactions are extracted. Conven- 
tionally, this is the FM stat o 10 ' 12 ' , but it is possible to employ also the AFM reference state or any other reference 
state such as the DLM state. In the DLM state the spin-spin correlations are missing while in the FM- or AFM-states 
the spins are strongly correlated. The DLM state is thus particularly suitable for discussion of the electronic origin of 
the magnetic phase stability as no specific magnetic order is assumed in this reference state. A traditional approach to 
the exchange interactions in the DLM-reference stat o n i 15 ' 16 rests on the generalized perturbation method (GPM)ji£ 
however, the final formula can also be obtained from the method of infinitesimal rotations of local moments applied 
to a random binary alloy. This equivalence is briefly sketched in Section II and in the Appendix. 

Once the exchange interactions are found, and it should be noted that they can be found even for very distant 
pairs of spins, one can study the ground state of the classical HH. In this way one can search over much broader class 
of possible magnetic systems, including incommensurate configurations, than it is possible by using a conventional 
supercell approach. It should be noted, however, that the supercell approach can be highly accurate and useful and 
we wish to point out the usefulness and complementarity of both above mentioned approaches for various aspects of 
the magnetic stability rather than to argue that one approach is better as compared to the other. 



The electronic properties of the Fe-overlayer on random Ta^Wi-^OOl) substrate will be studied in the framework of 
the tight-binding linear muffin-tin orbital (TB-LMTO) method combined with the coherent potential approximation 
(CPA) to treat alloy disorder in the substrate as well as the magnetic disorder in the DLM state. The effect of the 
surface is included in the framework of the surface Green function (SGF) approach^ which employs a realistic semi- 
infinite sample geometry (no slabs or periodic supcrcclls). The one-electron potential is treated within the atomic 
sphere approximation, but the dipole barrier due to the redistribution of electrons in the vicinity of the surface is 
included in the formalism. The TB-LMTO-SGF method can include the effect of layer relaxations approximately 
provided that they are known either from high-accuracy ab-initio calculations or from experiment. 

An important advantage of the TB-LMTO-SGF approach is the possibility to estimate exchange interactions be- 
tween magnetic atoms in the overlayer. We will describe such magnetic interactions in terms of the two-dimensional 
classical Heisenberg Hamiltonian, 



in which denotes the exchange integral between Fe-spins at sites i and j on the surface, and and e-,- are unit 
vectors in the directions of the local magnetization on sites i and j, respectively. 

Exchange interactions Jij can be evaluated in the the FM reference state and in the framework of the TB-LMTO- 
SGF method as 



II. FORMALISM 




(1) 




(2) 



where Aj(z) characterizes the exchange splitting of the Fe-atom at the site i, gfj(z) is the configurationally averaged 
Green function describing the motion of an electron between Fe-sites i and j in the magnetic overlayer on a random 
nonmagnetic substrate which corresponds to the spin a, a = (t,4-), and the integration is done over the contour 
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C in the complex energy plane z which starts below the valence band and ends at the Fermi energy. Symbol tr^ 
denotes the trace over the atomic orbitals (L = (l,m)). The quantity Aj(z) is defined in the TB-LMTO method in 
terms of the potential functions P°{z) of Fc atoms as Aj(z) = P}{z) — Pf{z). The above expression represents a 
straightforward generalization of the corresponding bulk expressio n 10 ' 13 since it is formulated in the real space. We 
refer the reader to Refs. [I8I - I20I for details of evaluation of the real-space configurationally averaged Green function 
9ij{z) for the magnetic overlayer on a random substrate. 

In the DLM-rcfcrcncc state, the exchange interactions Jij bear a standard form of effective pair interactions of the 
GPM. They are given explicitly by 

Jij = Im J c tTL { T *( z ) 9ij(z) Tj(z) gji{z)} dz , (3) 

where the gij(z) is the spin-independent configurationally averaged Green function in the DLM state between Fe-sites 
i and j. The electronic structure of the DLM state is treated in an alloy analogy, i.e., the Fe-sites are randomly 
occupied by two Fe-species with equal probability, representing local Fe-moments pointing in two opposite directions. 
The quantity r,(z) is defined as tj(z) = t\{z) — t\{z), where the symbols tf(z) (er =t,4-) refer to single-site t-matrices 
of the TB-LMTO method describing scattering due to both magnetic species on the i-th site with respect to a spin- 
independent effective DLM medium, see Eq. (|A.6|) . In the Appendix, the form of Eq. ^ is derived from the previous 
formula, Eq. ([2]), applied to the DLM state within the alloy analogy. 

The exchange integrals Jij characterize magnetic interactions between two particular atomic sites. It is also con- 
venient to define the exchange parameter J that reflects the molecular field experienced by a single moment in the 
ferromagnetic layer 

Jo = Yl J oi ■ ( 4 ) 

A negative value of this parameter indicates an instability of the FM state. A more detailed analysis of the magnetic 
stability employs the lattice Fourier transform of the pair exchange interactions Jij, see Section III C; the parameter 
Jo represents a special case of the Fourier transformed interaction, namely, its value for zero reciprocal vector, see 
Eq. ©. 



III. RESULTS AND DISCUSSION 



The calculations were done assuming the Vegard's law for the bulk lattice constant a of bcc-Ta^Wi-z random 
alloy (etTa = 3.300 A and a-yv = 3.165 A). The electronic structure is calculated selfconsistently for a model of the 
system which consists of 7 atomic layers representing the alloy substrate, one Fe layer, and 4 layers of empty spheres 
representing vacuum that are embedded between the semi-infinite substrate and the semi-infinite vacuum. Similarly as 
in Ref. [H, we assume a constant inward layer relaxation of 18 % between the Fe-overlayer and the bcc(001)-Ta 2: Wi_ 2 ; 
substrate. We have employed the spd/-basis and the Vosko-Wilk-Nusair exchange-correlation LSDA potential in all 
calculations. We used 20 points on the contour in the complex energy plane and 210 k-points in the irreducible wedge 
of the surface Brillouin zone in sclfconsistent calculations and up to 19600 k-points in the surface Brillouin zone to 
calculate exchange integrals for 172 inequivalent neighbors. 



A. Density of states 

We compare in Fig. [T] the Fe-projected densities of states (PDOS) for an unsupported monoatomic Fe(001) layer 
(assuming the in-plane lattice constant of bcc-W, but the results for that of bcc-Ta are very similar) with the PDOS 
for Fe/W(001) and Fe/Ta(001). Results for the FM and for the DLM states are shown. We observe a much stronger 
effect of hybridization among the Fe-located and the substrate states in the case of the W substrate than for the Ta 
substrate. The PDOS for Fe on the Ta substrate, especially that for the majority spin, is only weakly broadened as 
compared to the unsupported Fe layer in striking constrast to the case of the W substrate. We thus argue that the 
strong Fe-W hybridization is the reason for the AFM ordering on the W-substrate: the ground state of the weakly 
hybridized Fc/Ta(001) is similar to that of the unsupported Fe layer, namely the FM state. Roughly speaking, the 
stronger Fe-substrate hybridization also reduces the DOS at the Fermi energy in the non-magnetic state and thus 
reduces the tendency to fcrromagnctism (the Stoncr criterion) . It should be noted that the generally broader minority 
PDOS as compared to the majority one is due to the large exchange splitting of Fe electron levels: majority levels 
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are closer to the nucleus and thus are more strongly bound giving a narrower band and just the opposite holds for 
the more loosely bound minority bands. 

Another indication of a stronger hybridization of the Fe-W states as compared to the Fe-Ta states is documented 
by Fig. [3J where the concentration dependence of local magnetic Fe-moments evaluated in the FM and DLM states is 
shown. Generally, a stronger hybridization of Fe-W states results in an increasing difference between FM and DLM 
magnetic moments. The local magnetic Fe moments for the unsupported Fe layer are large and almost the same 
(about 3.35 /is) in both the FM and DLM configurations. A stronger hybridization reduces more the local magnetic 
moments for the Fe/W(001) overlayer as compared to the Fe/Ta(001) one. 

In Fig.|3]we present layer-resolved PDOS for first few top layers of the Fe/W(001) and Fe/Ta(001) systems (averaged 
over spins) and corresponding bulk substrates. We first note a large similarity between the bulk Ta and W total DOS 
which justifies the virtual crystal approximation used in Ref. [U to study the effect of alloy disorder in the substrate. 
The main difference is the position of the Fermi energy which is shifted downwards in the bcc-Ta because of a smaller 
number of valence electrons. As a consequence, the Fermi energy lies inside the bonding peak of the bcc-DOS while 
it is shifted to the energy region between the bonding and antibonding states in the bcc-W. We also mention that 
only the interface substrate layer and to some extent also the second substrate layer differ noticeably from the bulk 
substrate while the DOSs of other substrate layers are already quite similar to those of bulk. This is, of course, 
a typical feature of metallic systems in which the corresponding Friedel-like oscillations (which originate from the 
presence of an abrupt change of the charge density at the surface) are quickly damped inside the substrate due to the 
high density of screening electrons. 

B. Exchange integrals 

Concentration dependence of the exchange integrals J s for the first few shells of neighboring atoms s — 1 — 5 is 
shown in Fig. [4] First, this figure illustrates the principal difference between the approach used in Ref. [U and the 
present approach. In the former method, the total encgy differences between three magnetic configurations, namely 
the FM, c(2 x 2)- and p(2 x 1)-AFM configurations, are used to estimate two exchange integrals. One can call this 
approach the first-principle fitting as compared to the direct first-principle calculation of exchange integrals in our 
present approach. Clearly, realistic exchange integrals are not limited to the first two shells and, in particular, not in 
the two-dimensional case where exchange integrals can decay in general more slowly than in the three dimensional case 
(see e.g. Ref. [HH). It should be noted that we have plotted the exchange integrals without corresponding degeneracies 
(multiplicities of equivalent surface atoms in the shells) that are 4 or 8 in the present case. 

The most remarkable feature observed in Fig. [4] is a dramatic change of the leading first nearest-neighbor (NN) 
interaction from the strongly AFM coupling in Fe/W(001) to the FM coupling in Fe/Ta(001). The crossover between 
the AFM and FM couplings roughly coincides with the crossover between the FM and AFM ground state of the Fe 
overlayer. The concentration dependence of other interactions on the substrate alloy composition is generally weak. 
The second NN interaction has the AFM character in the whole concentration range while the third NN interaction 
is FM-likc for a W-rich substrate. The first two NN interactions agree reasonably well with those fitted from total 
energies for three magnetic configurations 2 - both qualitatively and quantitatively (note that the definitions of the 
Heisenbcrg Hamiltonian used in the present paper and in Ref. |2| differ by a factor of two) . The quantitative agreement 
of these NN interactions also proves that the chemical disorder in the TaW alloy substrate, which Ref. [U treated in a 
simple virtual crystal approximation but which the present study treats in the CPA, is really weak. 

In Fig. [5] we compare the difference between total energies of the DLM and FM phases with the calculated exchange 
parameter Jo, Eq. (|4]), for the DLM reference state adopted in this paper. Overall good agreement between both 
dependences over the whole concentration range should be pointed out. This confirms that the DLM state is a 
reasonable reference state for the estimate of exchange integrals. The small deviations are due to approximations 
inherent to the Heisenberg model (fixed size of local moments, the neglected contribution of induced moments in the 
substrate). 

The observed concentration trend of the first NN interaction and of the relative stability of the FM and DLM 
states can be understood qualitatively on the basis of general properties of a broad class of physical quantities in 
tight-binding models as functions of band filling] 21 ! 22 These quantities include, e.g., total-energy differences, local and 
non-local susceptibilities, magnetic and chemical interactions, etc.; all of them change their sign due to variation of 
the Fermi energy across the d-band. In the present case, the position of the Fermi level with respect to the Fe <i-band 
is directly controlled by the concentration of the Ta-W alloy substrate. 

Finally, in Fig. [5] we illustrate the convergence of the exchange parameter Jo, Eq. (U), as a function of the number 
of shells included in the summation. Note that now the shell degeneracies, N s , are included. Obviously, the inclusion 
of about 20 NN shells is sufficient for practical purposes (the deviation of the partial sum from the converged value 
is less than 1 %). At least 10 shells have to be considered in order to predict the formation of an incommensurate 
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magnetic ground state of Fe/Tao.2Wo.s(001) (see the next subsection for details). This is demonstrated in the inset 
of Fig. |51 where the lattice Fourier transform J(q) of the real-space exchange interactions Joi is magnified for the 
relevant wavevectors q in the surface Brillouin zone. The position of the minimum of the spin-spiral energy, i.e., the 
wavevector q for which J(q) has maximum, begins to approach the converged position for more than 10 shells while 
about 70 shells are needed to reach the converged position with the precision of 0.1 % of the wavevector q. The first 
two dominating NN interactions are obviously not enough to reproduce this local spin-spiral energy minimum. 

C. Magnetic stability 

Let us consider a class of possible arrangements (the so-called q— waves) of magnetic moments defined as 

e.; = (sin^ cos (/>i, sin^i sin</>j, cos^j) , (5) 

where the polar angle 9i = 8q is a constant and the azimuthal angle <pi = q • is a function of the position vector 
of site i. The energy per lattice site corresponding to the HH (QJ i?(q, #o) = — >^(q) sin 2 #o — J{0) cos 2 9o is expressed 
in terms of the lattice Fourier transform of the site-dependent exchange integrals 

•/•:q: X •/.•,<"' "• • (6) 

3 

It is easy to show that 

min.E(q, 6q) = — maxJ(q) (7) 
q,8o q 

and that 6q = ir/2 if the minimum of energy is achieved for q ^ 0. Let us note that J(q) is closely connected to 
the energies of magnetic excitations (magnons) in the system decribed by the HH, Eq. ([TJ. The maximum of J(q), 
or, equivalcntly, the minimum of — J(q), reached for a particular value of the vector q = qo in the surface Brillouin 
zone indicates a tendency of the magnetic system to form a magnetic ground state characterized by that wave vector 
qo- The wave vector qo = (point T in the surface Brillouin zone) correponds to the FM ground state while a 
nonzero wave vector qo corresponds to a more complex ground state including possible AFM or spin-spiral states. 
For example, qo at the point M of the surface Brillouin zone corresponds to the c(2 x 2) (checkerboard-like) AFM 
ground state and qo at the point X corresponds to the p(2 x 1) (row-wise) AFM ground state. A minimum of — J(q) 
outside any point of high-symmetry in the surface Brillouin zone suggests a tendency to form an incommensurate 
(spin-spiral) magnetic ground-state structure. 

In this way, we are able to investigate the stability of a much broader class of magnetic systems as compared to 
the conventional total energy search. It should be noted, however, that possible more general ground states are not 
included in the present search and a more advanced search based on the ovcrlayer Heisenbcrg Hamiltonian, Eq. ([lj, 
would be needed. Finally, we mention that in the present study we neglect the anisotropic part of the effective 
magnetic Hamiltonian, i.e., the anisotropy that concerns orientation of the spins with respect to the underlying 
lattice. Thus the relativistic effects (the magnetic anisotropy and the Dzyaloshinskii-Moriya interactions (DMI)) as 
well as the magnetostatic (magnetic dipole-dipolc) interactions are missing in the present analysis. The dipole-dipolc 
pair interactions between the first and second NN Fe-Fe pairs are at least two orders of magnitude weaker (below 0.01 
meV) than the typical exchange interactions in the present case. However, their long range contributes significantly 
to formation of domain walls and it might influence other details of the magnetic order too. It is also known that the 
DMI can lead to long-range chiral structures due to weakening of FM-interactions, as it happens in the Mn/W^lO) 2 ^ 
and Mn/W(001)2i systems. We note that in a consistent theoretical treatment all possible anisotropic terms would 
have to be included on the same footing. 

The curves of — J(q) plotted in Fig. [7] and especially the positions of their minima reflect the dependence of the 
magnetic ground state on the composition of the substrate. A clearly pronounced minimum for Fe/W(001) at the M 
point in the surface Brillouin zone confirms the total energy search, namely that the checkerboard c(2 x 2)-AFM is 
the ground state of the system. With increasing Ta content, the stability of the c(2 x 2)-AFM state weakens and for 
about 20 % of Ta atoms a weak minimum corresponding to an incommensurate state (with the wavevector on the line 
X — M close to M) develops. It should be noted that this minimum is sensitive to the number of shells included in the 
evaluation of — J(q) and is missing if we consider only few leading exchange interactions (see previous subsection). We 
cannot exclude the possibility, however, that this shallow local minimum will disappear when the relativistic effects are 
included. The magnetocrystalline anisotropy, for example, strengthens the tendency towards collinear AFM alignment 
and competes with the DMI interaction which promotes the rotating, non-collinear magnetism as it was demonstrated 
recently for the Fe-double chains on the fcc-Ir(OOl) surfaced. If the Ta content further increases, the p(2 x 1)-AFM 
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configuration becomes the ground state as it is illustrated for the case with 50 % of Ta atoms where the minimum 
of — J(q) occurs at X, see Fig. [7] In the Ta-rich alloy the stability of the FM phase increases as compared to the 
c(2 x 2) one until for Fe/Ta(001) the wave vector at T (or the FM state) becomes the magnetic ground state. 

The observed tendency to a non-collinear ground state for x « 0.2 can be explained by frustration effects on the 
square Fe lattice accompanying the dominating first and second nearest-neighbor interactions that are both negative 
for concentrations x < 0.45, see Fig.|U For W-rich substrates (x — > 0), the first nearest-neighbor interaction dominates 
(| Ji | 3> | J2I) and the c(2 x 2) AFM ground state is not frustrated. For equiatomic concentrations (x — > 0.5), the second 
nearest-neighbor interaction is the strongest one (IJ2I ^ \ Ji\) and the p(2 x 1) AFM ground state is not frustrated 
cither. However, for compositions around x ~ 0.25, both interactions are of comparable magnitudes ( J\ « J 2) which 
leads to pronounced frustration of these simple AFM states. Moreover, these AFM states become nearly degenerate 
at x = 0.2, where the two leading interactions satisfy roughly a relation Jj = 2J 2 £ The formation of a single-q spin 
spiral for x ~ 0.2 (or of a more complex spin arrangement not considered here) represents thus a natural consequence 
of existing interactions in the Fe-monolayer. 

Finally, the present results (e.g., the critical concentration for the AFM to FM crossover) can be influenced quan- 
titatively by possible structural changes in real system. A recent study^ indicates a possibility for the B2-ordering 
in a bulk bcc-TaW alloy over a broad concentration range. Also, a smaller surface energy of bcc-Ta as compared to 
the bcc-W— together with the larger size of Ta atoms seems to indicate a possibility for the segregation of Ta-atoms 
at the bcc-(Ta,W)(001) alloy surface. On the other hand, the qualitative conclusions of the present paper remain 
unchanged. 



We have developed a new approach to study the magnetic phase stability of magnetic overlayers on non-magnetic 
substrates. The approach consists in the evaluation of exchange integrals between local magnetic moments in 
the magnetic overlayer using the adiabatic approximation and the real-space Green function approach. This 
allows us to use this approach also for magnetic overlayers on a random substrate and/or for random alloy 
magnetic overlayers. Estimated exchange interactions between pairs of local moments in the overlayer are used 
to construct the effective two-dimensional Heisenberg Hamiltonian, whose stability with respect to periodic spin 
excitations is investigated. The maxima of the lattice Fourier transform of the site-dependent exchange integrals 
are searched for as indications of stable periodic spin structures. Such a search allows one to investigate the 
stability of a much broader class of possible magnetic configurations as compared to the conventional total energy 
search which is limited to a few empirically chosen configurations. The present approach can be extended in 
future to include more general magnetic configurations and/or anisotropic effects (spin-orbit and dipole-dipolc 
interactions); the latter are indispensable for estimation of the corresponding Curic/Neel temperatures. As a case 
study we have investigated in detail the magnetic phase stability of the Fe overlayer on the bcc-Ta x Wi_ 2; (001) 
random substrate. The exchange interactions for this system were extracted from the DLM state, which involves 
no correlation among spins of overlayer atoms. Our results are in good agreement with the results of a recent 
study based on the total energy search and confirm basic approximations adopted in this study. In addition, we 
have predicted a possible incommensurate magnetic configuration for the W-rich substrate alloy (at about 20 % of Ta) . 
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The purpose of this Appendix is to show a close relation of the formula ([3]) for exchange interactions in the DLM- 
reference state to the well-known formula ([5]) valid for any collinear reference state, including the FM state. The 
latter formula for the pair interaction between local moments on lattice sites i and j can be written as 



IV. CONCLUSIONS 
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Appendix: Method of infinitesimal rotations in the DLM state 




(A.l) 
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where the integrated function in the TB-LMTO method has a form 



Wii = tr L 



{(?! 
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(A.2) 



where the energy argument z of each quantity has been omitted. Here we assume for simplicity that possible chemical 
disorder (random occupation by different species) is confined to non-magnetic sites whereas the sites i and j in ([A. 11 
IA.2I) arc occupied by a single species with non-zero local moment. 

The standard treatment of the DLM state employs the CPA2£ which leads to the alloy analogy, i.e., the electronic 
structure can be obtained from a binary equiconcentration random alloy A0.5B0.5 with collincar spin structure, where 
the spin-dependent potential functions of both components are given by P^ a = Pf-° = P? (where a =|, |, and 
— 1=4-, — 4-=t)- The symbol A (B) represents atoms with local moments pointing upwards (downwards); for symmetry 
reasons, the configurationally averaged Green function as well as the coherent potential functions are spin-independent, 
i.e., gfj = g i:j and T? = V l . 

Exchange interactions in the DLM-reference state can naturally be obtained from energy changes accompanying 
infinitesimal rotations of local moments that were originally pointing in the same direction. This corresponds to 
occupation of both sites i and j by the same atomic species, e.g., by A atoms. The configuration average of the 
exchange interaction leads to its original form (|A.1[) with the function uiij(z) replaced by 



tr 



{( 



R 



At 



P 



Al\ =A,Af 



pAt _ pM 



-A.AJ. \ 
Jji J ) 



(A.3) 



where g^' a denotes the conditionally averaged Green function in spin channel a between sites i and j occupied by 
atoms A. Note that the conditionally averaged Green functions depend on the spin index. Their values within the 
CPA and for i ^ j are given by 



-A,A<r 
iij 



(A.4) 



where the spin-dependent factors on the r.h.s. can be expressed in terms of the spin-independent on-site Green 
functions ga and coherent potential functions Vi as 

rAa r-i 1 / 7~>cr t-> \ - 1 — 1 
Ji 
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[l + (P?-Vi)gu} 
[1 + 9u (P? - Vi)] 
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The scattering due to the atomic species A and B with respect to the effective DLM-medium is described by 
spin-dependent single-site t-matrices tf a = tf~ a = tf , where 



% = ft{P?-'Pi) = {P?-Vi)ft. 
Their difference can be related to the difference of the potential functions Pf, namely, 
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and, similarly, 



(A.7) 



(A. 



Substitution of Eq. (|A.4I) into Eq. (|A.3[) . cyclic invariance of the trace and identities (|A.71 IA.8[) lead to the final 
expression of the function Wij(z) in the pair interaction Jij (IA.1I) : 



= tr L 



(A.9) 



This result is equivalent to the formula ^ in the main text. 
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FIG. 1: (color online). Spin-resolved densities of states projected on Fe-atoms for various systems in the ferromagnetic and 
DLM state: (a) an unsupported Fe monolayer, (b) Fe monolayer on the top of Ta(OOl) surface, and (c) Fe monolayer on the top 
of the W(001) surface. The majority (minority) spin states are shown in the upper (lower) panels. The vertical lines denote 
positions of the Fermi levels. 
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FIG. 2: (color online). The local magnetic moments mp e of Fe atoms in the monolayer on Ta^Wi-^OOl) substrate as a function 
of the alloy composition: (a) the DLM configuration (crosses, |rriFe| is shown) and (b) the FM configuration (circles). 
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FIG. 3: (color online). Layer-resolved densites of states averaged over spins for the Fe-overlayer in the disordered local moment 
state (dotted-dashed line), first four substrate layers (full line), and for the bulk substrate (dashed line): (a) Fe/Ta(001) system, 
and (b) Fe/W(001) system. The curves are shifted vertically with respect to each other. Note the different scale for (a) and 
(b) plots. The vertical lines denote positions of the bulk Fermi levels. 
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FIG. 4: (color online). Exchange integrals between Fe-atoms (up to the fifth shell) in the Fe overlayer on the Ta^Wi-^ (001) 
random alloy surface as a function of the substrate alloy composition. 



60 




-60 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

W <- y -> Ta 



FIG. 5: (color online), (a) The total energy difference of the disordered local moment (DLM) and ferromagnetic (FM) 
configurations per surface atom (-E°t M — ^tot > solid circles). Negative values indicate a tendency to the antiparallel alignment 
while the positive values indicate a tendency to the parallel alignment of local magnetic moments, (b) Exchange parameter Jo 
in the reference DLM configuration defined as the sum of corresponding exchange interactions Joi in the Fe-overlayer (crosses). 
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FIG. 6: (color online). The partial sum of the exchange integrals, Jo(smax) = X]J™i x N S J B for the Fe/Tao.2Wo.8(001) system 
as a function of the shell number s max - Js is the exchange integral for the atomic shell s and N a is the corresponding shell 
degeneracy, i.e., the number of equivalent atoms in a given shell. The dependence of the lattice Fourier transform — J(q) (see 
Fig. Efb)) on the number of shells s max used in the Fourier transform is shown in the inset. The last of the curves in the inset 
(s ma x=100) is essentially a converged result with respect to the shell number. 




FIG. 7: (color online). The lattice Fourier transform J(q) of exchange integrals in the Fe overlayer as defined by Eq. (jGj along 
a chosen path in the first surface Brillouin zone plotted for different substrate alloy compositions: (a) bcc Fe/W(001), (b) bcc 
Fe/Tao.2Wo.8(001), (c) bcc Fe/Tao.5Wo.s(001), and (d) bcc Fe/Ta(001). The minimum of — J(q) for each substrate composition 
is indicated by a symbol on the corresponding curve. The inset shows the first surface Brillouin zone and the corresponding 
high-symmetry points. 



